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Abstract. We examine the plausibility of crossing the cosmological constant (A) barrier 
in a two-field quintessence model of dark energy, involving a kinetic interaction between the 
individual fields. Such a kinetic interaction may have its origin in the four dimensional effective 
two-field version of the Dirac-Born-Infeld action, that describes the motion of a D3-brane in 
a higher dimensional space-time. We show that this interaction term could indeed enable the 
dark energy equation of state parameter to cross the A-barrier (i.e., = —1), keeping 
the Hamiltonian well behaved (bounded from below), as well as satisfying the condition of 
stability of cosmological density perturbations, i.e., the positivity of the squares of the sound 
speeds corresponding to the adiabatic and entropy modes. The model is found to fit well with 
the latest Supernova Union data and the WMAP results. The best fit curve for w,^ crosses 
— 1 at red-shift z in the range ~ 0.215 — 0.245, whereas the transition from deceleration to 
acceleration takes place in the range of z ^ 0.56 — 0.6. The scalar potential reconstructed 
using the best fit model parameters is found to vary smoothly with time, while the dark energy 
density nearly follows the matter density at early epochs, becomes dominant in recent past, 
and slowly increases thereafter without giving rise to singularities in finite future. 
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1. Introduction 

A variety of recent observational probes, including in particular the type la Supernovae (SN la) 
[1, 2, 3, 4, 5, 6, 7, 8, 9], indicate that our universe has entered in a phase of accelerated expansion 
in recent past, following an early decelerating regime. Despite several alternative proposals, 
such as modified gravity [10] and the averaging of cosmological inhomogeneities [11], the origin 
of this acceleration has widely been attributed to a 'mysterious' energy component, namely 
the dark energy (DE), which constitutes about 72% of the present universe. Moreover, the 
cosmic microwave background (CMB) temperature fluctuation measurements by the Wilkinson 
Microwave Anisotropy Probe (WMAP) [12, 13, 14] as well as the large scale red-shift data from 
the Sloan Digital Sky Survey (SDSS) [15] indicate that our universe is very nearly spatially 
fiat, so that spatial inhomogeneities may be neglected at large scales. Although the DE 
closely resembles a positive cosmological constant A, for which the DE equation of state (EoS) 
parameter — —1, there are some serious theoretical problems, such as fine tuning and 
coincidence, associated with A [16] . Specifically, if the DE is supposed to be due to A and the 
acceleration began only in recent past, then (i) what makes the DE density scale very small 
compared to the Planck scale? and (ii) why is the DE density is of the order of the present 
critical density P(,^ right now? Hence, there have been suggestions that the DE may be (more 
appropriately) dynamic and can be modeled by one or more scalar ficld(s) originating from 
a fundamental theory. Of major interest arc the DE models developed in the framework of 
quintessence and tracker fields [17, 18], k-essence [19, 20], Chaplygin gas [21] etc., (see [22] for 
extensive reviews). However, in many of these models the value of is always restricted to be 
> — 1, which is not desirable for a consistent statistical fit with the observational data. In fact, 
even with the presumption that is a constant, the recent WMAP five year data, combined 
with with those for SN la and baryon acoustic oscillation (BAO) peaks, constrain the value of 
1 +w^, to be between —0.14 and 0.12, at 95% confidence level (CL) [14]. For a time-varying 
DE, the same data constrain the value w^^^ of the DE EoS parameter at the present epoch 
(i.e., at red-shift 2; = 0) to be between —1.33 and 0.79 (at 95% CL) [14]. Since in the distant 
past, the value of a variable must have to be ^ —1 (so that the universe had a decelerated 
expansion and structures were formed), there is a fair plausibility that one (or possibly more) 
transition (s) from > — 1 to < —1 (or vice versa) could have been taken place in the 
recent course of evolution of the DE, and at present w-^ — w^^ < — 1. 

The crossing of the cosmological constant barrier {w-^ — —1) can, most simply, be achieved 
in the so-called quintom scenario [23], where there are two (or more) scalar fields, (at least) 
one of which is of 'phantom' nature, i.e., carries a wrong sign in front of the kinetic term in 
the Lagrangian [24]. Such a phantom field is quantum mechanically unstable [25] and also 
gives rise to singularities in finite future [26, 27]. Moreover, classical instabilities could arise as 
the dominant energy condition gets violated in the models involving the phantom fields [28]. 
Attempts have therefore been made to circumvent the problem of = — 1 crossing in various 
alternative ways. Notable among these are the scalar-tensor models [29], brane- world models 
[30], multi-field k-essence models [31, 32], modified gravity models [33], string-inspired dilatonic 
ghost condensate models [34], quantum-corrected Klein-Gordon models with quartic potential 
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[35], coupled DE models [36], H-essence (complex scalar) models [37], etc. However, apart from 
a very few exceptions (such as the scalar-tensor models [29], or models where the kinetic term 
abruptly flips sign due to some extraordinary nature of the potential [38]) the = — 1 crossing 
is hard to be realized with a single-field DE. Even in the case of a single-field k-essence DE, 
with a generic non-linear dependence of the Lagrangian on the kinetic term, such a crossing 
either leads to instabilities against cosmological perturbations or is reahzed by a discrete set of 
phase space trajectories^ [40]. In multi-field DE models, however, the — —1 crossing could 
be made possible, as is shown for example in refs. [31, 32], although the field configuration 
may be severely constrained by the criterion of stability, i.e., the square of the effective speed 
of propagation of cosmological perturbations should be positive definite [32]. 

In this paper we explore the plausibility of the A-barrier crossing in the framework of 
a two-field quintessence model with a kinetic interaction between the individual fields. Such 
a model may be looked upon as a specialization of a more general (interacting) multi-field 
k-essence scenario, which involves non-canonical (higher order) kinetic terms for the scalar 
fields [19, 32, 41, 42, 43]. Moreover, the kinetically interacting double quintessence (KIDQ) 
Lagrangian may, under certain approximations, be derived from the four dimensional effective 
two-field version of the Dirac-Born-Infeld (DBI) Lagrangian describing the evolution of D3- 
branes in higher dimensional string theoretic manifolds [44] . The biggest advantage with such 
a Lagrangian, compared to those in other A-barrier crossing multiple k-essence models [31, 32], 
is that the total DE Hamiltonian consists of a positive definite kinetic part, which ensures that 
it is bounded from below and the model is quantum mechanically consistent. Stability against 
cosmological density perturbations further requires the squares of the effective (sound) speeds 
of propagation of the adiabatic and entropy modes to be positive definite as well. For the DBI 
multiple scalar fields in homogeneous cosmological backgrounds, both these sound speeds turn 
out to be the same, implying isotropic propagation of the adiabatic and entropy modes [41, 42]. 
Assuming this result to hold approximately for KIDQ (which is an approximation to the DBI 
two-scalar scenario), we find the square of the effective (isotropic) sound speed to be positive 
definite, ensuring the stability of the KIDQ model§. 

We consider certain specific ansatzc to solve for the KIDQ field equations, and obtain the 
condition under which the = — 1 line could be crossed in some regime. In choosing the 
ansatze, we particularly emphasize on the following: 

(i) the kinetic energy densities of the interacting scalar fields should always be positive definite, 

(ii) the DE density should be less but not very smaller than the matter density at early epochs, 
and should dominate the latter at late times, and 

(iii) the DE density should not grow rapidly with increasing scale factor a (i.e., decreasing 
red-shift z) and reach to abnormally high values in finite future. 

X Of course, there arc exceptions as well, sec for example ref. [39]. 

§ More precisely, however, there is a splitting between the propagation speeds of the adiabatic and the entropy 
modes, when the KIDQ is taken to be an exact theory (not an approximation to DBI). This we find in a 
subsequent paper [45] (in preparation) by carrying out the stability analysis for KIDQ, following the general 
formalism worked out in refs. [42, 43] in the context of multi- field DBI and k- inflation. The squares of the 
propagation speeds turn out to be positive definite anyway. 
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These are important in order to avoid (i) ghosts or phantoms, (ii) coincidence or fine-tuning 
related problems, and (iii) occurance of future singularities, respectively. 

We then constrain the parameters of the model with the latest Supernova la data compiled 
in ref. [8], viz., the 307 Union data-set, as well as with the WMAP 5-year [14] update of the 
CMB-shift parameter TZ and the scalar spectral index n^, which determines the BAO peak 
distance parameter A from the SDSS luminous red galactic distribution [15]. After uniformly 
marginalizing over the Hubble constant Hq, we obtain good fits of the model with the data. 
The minimized value of the total (SN-I-CMB-I-BAO) is found to be ~ 311, which is better 
than the minimized x^(— 313) found with the Union data-set in ref. [9] for the cosmological 
constant DE coupled with cold dark matter - the so-called ACDM model. The best fit values 
of the parameters of our KIDQ model indicate that the crossing from > — 1 to < — 1 
takes place at a red-shift range 0.215 < Zc< 0.245, whereas the transition from the decelerated 
regime to the accelerated regime occurs in the range 0.562 < Zt < 0.603. At the present epoch 
{z = 0), the best fit values of the matter density parameter and the DE EoS parameter, are 
respectively found to he within the ranges 0.279 < f2om < 0.281 and —1.123 < w^^ < —1.077. 
All these results are fairly in agreement with those found with other model-independent or 
model-dependent parameterizations of the DE in the literature [46, 47, 48, 49, 50]. 

Finally, we integrate the scalar field equations of motion and reconstruct the interacting 
double quintessence potential using the best fit model parameters. We show that the 
reconstructed potential has a smooth dependence (i.e., without any discontinuity or multi- 
valuedness) on the scale factor a. We work out the approximate analytic expressions for the 
potential as function of the scalar fields, and find that they also exhibit the same smooth nature 
at early and late stages of the evolution of the universe. 

This paper is organized as follows: in sec. 2 we describe the general framework of the multi- 
scalar (k-essence) DE scenario, following the formalism shown in refs. [42, 43] in the context of 
multi-field k-infiation. In sec. 3 we emphasize on a special case which involves two quintessence 
type of scalar fields with canonical kinetic terms in the Lagrangian, and with a specific kinetic 
interaction between the individual fields. Assuming suitable ansatze for the solutions of the field 
equations we work out the condition under which the cosmological constant barrier = — 1 
could be crossed, and find the expression for the Hubble parameter maintaining this condition. 
In sec. 4 we fit our KIDQ model with the 307 Union SN la data [8], combined with the 
CMBH-BAO results from WMAP and SDSS, to obtain the DE density and EoS profiles. In 
sec. 5 we use the best fit values of the model parameters to reconstruct phenomenologically the 
interacting double quintessence potential and determine the temporal variations of the scalar 
fields. We also work out the approximate analytic functional forms of the potential in terms of 
the scalar fields, at early and late stages of the evolution of the universe. In sec. 6, we conclude 
with a summary and some open questions. In the Appendix, we show how the KIDQ action, 
that we consider, could be derived from the two-field DBI action under certain approximations. 
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2. General Formalism 

Let us consider the following action, in (3 + 1) dimensions, for gravity minimally coupled with 
matter fields and N number of kinetically interacting (k-essence) scalar fields (p^ {I — 1, . . . , N): 

d'xV^l^ + Cm + P{X''A'')] , (1) 



2k 

where — SttG is the gravitational coupling constant, is the Lagrangian density for matter, 
that is considered to be pressureless dust. P{X^^ , (f)^) is the multi-scalar Lagrangian density, 
with 

X'' ^ - \ g^' d,<t>' dA\ (/,J=1,...,7V), (2) 

describing the kinetics of the scalar fields [42] . 

In a spatially fiat Priedmann-Robertson- Walker (FRW) background, with fine element 

ds" = -dt^ + (t) [dr^ + r^ [d^^ + sin^ M^^) ] (3) 

the above expression for X^"^ reduces to 

x^' = x'^ - ^<P''<P'', (4) 

where the dot denotes time derivative (d/dt) and the prime denotes derivative (d/da) with 
respect to the scale factor a, which has been normalized to unity at the present epoch t — tg. 
H = d/a is the Hubble parameter. 

The Friedmann equations and the scalar field equations of motion are given by 

f!) =t(^-+^-^ ' H=^-^^-'^[p^ + {p^+p^)] , (5) 
\a J 6 a 1 

d { . dP -A ^ dP 

where is the energy density of matter in the form pressureless dust, and Px-iPx 
multi-field dark energy density and pressure, given respectively as 

BP 

Px = 2X^'^ - P ^ P. - P- (7) 
Assuming that there is no mutual interaction between matter and dark energy, the 
Friedmann equations (5) integrate to give p^ — Po^a~^, where p^^ is the matter density at 
the present epoch {t — to,a — 1). One also has the continuity equation for the dark energy 

Px^-3Hip^+p^) p'^ = _^(p^+p^) . (8) 

From the Friedmann equations (5) one obtains the expressions for the DE EoS parameter 
the total EoS parameter w, and the deceleration parameter q: 
Px . 2X^-^ dP 

d l + 3w , , 
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where p^^ = 3Hq/k^ is the present critical density; Hq being the value of H at the present 
epoch {t — to). 

= + (12) 

is the normalized Hubble parameter and flom — Pom/Poc is the present matter density parameter. 

The transition from the decelerating regime to the accelerating regime takes place when the 
deceleration parameter q changes sign, i.e., the total EoS parameter w becomes less than —1/3, 
by Eq. (11), and the DE EoS parameter is further less, by Eq. (10). The crossing from 
W-^ > — 1 to < —1, on the other hand, requires a flip of sign of the quantity X^'^dP/dX^'^ , 
presuming that the DE density is positive definite. In the next section, we examine the 
plausibility of such a crossing by considering for simplicity a model involving only two fields 
{N = 2) with usual (canonical) kinetic terms (quintessence type), but with a specific type of 
kinetic interaction, which could have its origin in the two-field DBI action, as we show in the 
Appendix. 



3. Kinetically interacting double quintessence 

Let us take into account the following special form of the Lagrangian density for the DE, 
consisting of only two scalar fields: 



P = 5uX'' - - f (<^/,J-i + 5j.,^j)X^J - y(0O, (13) 

where 7 are positive constants, V{(j)^) is the scalar potential, and the indices /, J run for 1, 2. 
Denoting the two fields as 0^ = (0,0) we can re- write the above Lagrangian as 



^ = Y + f - iQiM - V{4>,^), where Q{4>A) = yl - ^ 4>t (14) 

This implies that the scalar fields and £ have usual (canonical) kinetic energy densities 
(given respectively by the first two terms on the right hand side), and therefore are similar 
to ordinary quintessence fields. However they have a mutual kinetic interaction of a specific 
form proportional to (^(0, given above, which may originate from the two-scalar DBI action, 
approximated for /3 ^ 1 and 7^1 (but 7"^ <S /3) as shown in the Appendix. 

The dark energy pressure is equal to P in Eq. (14), whereas the expression (7) for the 
dark energy density reduces to 

The DE equation of state parameter w^, Eq. (9), now takes the form 



^ = - 1 H 

Px Px 



+ e + 



(16) 



2Q(0,O. 

The presumption that the parameter /3 <^ 1, is in support of the positivity of the term 
under the square root in the expression for (^(0, given in Eq. (14). That is, the requirement 
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Q'^{4>,C) > for the vahdity of the model, could be fulfilled when /3 -C 1, even if and ^ 
vary fairly rapidly with time and the product (pC, > 0. Considering further, Q{(t>.^) itself to be 
positive, the kinematical part of the DE density, given by the first three terms (kinetic energy 
densities of the fields plus their kinetic interaction) on the right hand side of Eq. (15), remains 
positive definite. As such, the total DE Hamiltonian is bounded from below and the model is 
quantum mechanically consistent. Moreover, since /3,7 and Q{4^,C) are all positive, it follows 
from Eq. (16) that, < —1 (in some regime) necessarily implies the product 0^ < 0. In other 
words, the condition for the crossing of the — —1 barrier at a particular epoch, is that one 
of the two fields (0, ^) must fall off with time, whereas the other one should increase with time. 

Now, using Eqs. (15), (16) and the continuity equation (8), one obtains the following 
expression for the potential y as a function of the scale factor a: 



V{a) = 



7 



Q{a) 



+ A 



3 

a 



4>\d)+e{a) + 



2 Q{d) 



,(17) 



where A is an integration constant. Plugging Eq. (17) back in Eq. (15) we get the DE density 

Pi 0(a) i{d) 



as a function of a: 



A 



31 ^ 



d)+e{~a) + 



(18) 



2 Q{a) 

One may note that Eq. (17) could also have been obtained by using the scalar field equations 
of motion (6), which in the present scenario reduce to 



d 
~dt 

d 
dt 



+ 



^ + 



4g(0,e), 

/37 



4g(0,o, 

Under a dimensional re-scaling: 
. e 







dV 



dV 



A ^ 



the DE density, pressure, and the scalar potential change as 



7 

Poc 



Px Px 
Px^—, Px^ — 



Poc 



(19) 



(20) 



(21) 



Poc Poc 

while all the above equations (15) - (19) remain invariant. On the other hand, the expression 
(12) for the normalized Hubble parameter reduces to 



Om 



A + 



Om 



dd 

- 3 ^ 



Pi 0(a) i{a) 
2 Q{d) 



(22) 



Let us now consider the following ansatze for the kinetic energy densities of the scalar 
fields: 



(a) = - 0^(a) 



Pf(a) 



/(a) + V/'(a) - F 
/(a) - V/2(a) - fc2 



(23) 
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where /(a) is taken to be a positive definite and well-behaved function of a, A; is a positive 
constant, and /(a) > k at all epochs. Eqs. (23) imply that 

02 + = 2/(a), and 4> t ^ ±k. (24) 

We choose to take 4>C = so that the DE EoS parameter w^, Eq. (16), could be made less 
than —1 in some regime. Moreover, this choice guarantees the positivity of the square of the 
kinetic interaction, which now reduces to a constant: 



= 1 + f • (25) 
The expressions for the time derivatives of the scalar fields are given by 

v2 v2 
whereas from Eqs. (16), (18) and (17), we respectively obtain the following expressions for the 
DE EoS parameter and density, and the scalar potential: 



w^{a) — — 1 + 



- -2Q 



(27) 



pM = A + ^ Ina - ey -^da, (28) 

Via) ^A+'-^lna-m-l-eflfd-a. (29) 

Let us now assume a specific form of the function /(a), given by 

/(a) = Aa-" + k, where A>0, < u < 3 , (30) 

so that the criterion f{a)>k>0 is automatically satisfied. Furthermore, < < 3 ensures 
that /(a), and hence the kinetic energy densities and of the scalar fields, fall off with 
increasing values of the scale factor a. However, these fall offs are not faster than that of the 
matter density (p^ ~ 1/a^). This is essential in order that the quantities and 1^^, which 
compose the total DE density p^, come to dominate at late times, i.e., for large values of a. 
Eqs. (26) reduce to 

= > i\P) = , (31) 

and the Eqs. (27) - (29), for w^^p^ and V ^ take the form 

+ 6filna + A, (33) 





= - 1 + 


Pxia) 




= — a 


u 


V(a) 





0-1^ A(a-''-l)+ 6Slna, 



(34) 



where we have defined 



B — k [ — — 1 ] = constant , (35) 
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and Vq is the value of the scalar potential V at the present epoch {t — tQ,a — l): 

Prom Eqs. (22) and (33) one also obtains the following expression for the normalized Hubble 
rate H = H/Hq. 



H^a) ^ M + ^ + 651na + A. (37) 



At a = 1 (present epoch), H — 1, whence 

A = 1 - Qom - — , (38) 
and the above expression (37) reduces to 



HHa) = 1 + — (l-an+ ^(l-a')+ 6Slna. (39) 

In the next section, we fit this Eq. (39) with the latest Supernova la data [8], as well as 
with the CMB+BAO results from WMAP and SDSS [14, 15], and determine the DE density 
and EoS profiles over the red-shift range that is probed. 

4. Observational constraints 

Wc perform a analysis so as to constrain the model parameters A, B and r2om.) '^oi two specific 
choices of the index u (= 1, 2) in the ansatze (30). The SN la Union data-set [8], which we use, 
consists of 307 most reliable data points that range up to red-shift z — (1/a — 1) ~ 1.7, and 
include large samples of SN la from older data-sets [1, 2, 3, 4], high- 2; Hubble Space Telescope 
(HST) observations and the SN Legacy Survey (SNLS) [5]. 

The SN la data provide the observed distance modulus /jtobsi^i): with the respective la- 
uncertainty ai{zi), for SN la located at various red-shifts Zi,{i — 1, . . . , 307). The for the 
SN observations is, on the other hand, expressed as 

XsAt^oAkm.A.B) - --2T-T , (40) 

where 

li{zi) = 51ogio [L'i(zi)] +/io, (41) 
is the theoretical distance modulus. 

Dl{z.) - (1 + -.) r , , (42) 

is the Hubble free luminosity distance in terms of the parameters {flom, A, B), and 



Ho = 5 log 



10 



Mpc 



+ 25 = 42.38 - 51ogio/i, (43) 



h being the Hubble constant Hq in units of 100 Km s ^ Mpc ^. The parameter fiQ is a nuisance 
parameter, independent of the data points, and has to be uniformly marginalized over (i.e.. 
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integrated out). For such a marginalization one may follow the procedure shown in refs. 
[32, 51, 52], where is first expanded suitably in terms of /iq. Then one finds the value 
of /xq for which such an expanded form of is minimum. Substituting this value of /iq back 
in X^,j^, finally enables one to perform the minimization of the resulting expression with respect 
to the parameters iytom-i A, B), in order to determine the values of the latter best fit with the 
SN la observations. 

The CMB shift parameter TZ, that relates the angular diameter distance to the last 
scattering surface (at red-shift zis) with the co-moving sound horizon scale at recombination 
and the angular scale of the first acoustic peak in the CMB temperature fluctuations power 
spectrum [15, 53], is given by 

n{z,) = Q}£ r__^! , (44) 

Jq H{z; Qorn, ^, B) 

where z^, is the red-shift of recombination. The WMAP flve year data [14] updates z-^ — 
1090.04 ±0.93 and the observed shift parameter TZobsiz*) = 1.710 ±0.019. The for the CMB 
observations is given by 



XL» = ^5^^^, (45) 



where ctr is the la error in the WMAP data [14]. 

Now, the scalar spectral index n^, which determines the observed value of the BAO peak 
distance parameter Aobs from the distribution of the SDSS luminous red galaxies [15] through 
the relation Aohs = 0.469 (n,/0.98)"°-^^ ± 0.017, is updated by the WMAP five year data as 
Us = 0.960 ± 0.013 (see the first ref. of [14], see also ref. [52]). The theoretical expression for 
the distance parameter is, on the other hand, given by 

-1 2/3 



A ^ n, 



1/2 
Om 



1 r" dz 



(46) 



Zh^Hizh) H{z;flom,AB) 
where = 0.35. The x^ for the BAO observations is expressed as 

2 _ {Aobs - -A) , , 



a A being the la error in the SDSS data [15]. 

The total x^; which needs to be minimized in order to determine the likelihood of the 
model parameters (Qom, A, B) with the entire SN+CMB+BAO data, is thus given as 

Xtotal ~ Xgjv + Xc.jvfS ~^ ^BAO ' 

Of course, x^jy has already been minimized with respect to the nuisance parameter Hq, Eq. 
(43), by the process discussed above. 

For two specific choices of the index i> (= 1,2), that appears in the ansatz (30), the best 
fit values the parameters {flom, A,B), as well as the minimized value of Xtotah shown in the 
table 1. Fig. 1 shows the evolution of w^{z) and p^i^) (alongwith the corresponding la errors) 
throughout the entire red-shift range < 2; < 1.75 of the available data, for both the choices 
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Index 


Best fit 


model parameters 


Minimized 






A 


B 


Xtotal 


1 


0.2790 


0.2062 


0.2506 


311.07 


2 


0.2816 


0.0505 


0.0782 


311.27 



Table 1. Values of the parameters (Oo^^^j-B) of the model, best fit with 
SN+CMB+BAO observations, and the minimized total x^, for the choices u = 1,2. 




Figure 1. Evolution of w^{z) and p^{z) (best fit with the SN+CMB+BAO 
observations) throughout the red-shift range < z < 1.75) alongwith the corresponding 
Icr error (shaded) regions, are shown for the choices v = 1 (upper panels) and v = 2 

(lower panels). The point Zc denotes the red-shift at which the W-^ = — 1 line is crossed, 
and , are respectively the values of , at the present epoch {z = 0). 

of u. The maximum likelihood of the present value w,,^ of the DE EoS parameter is found to 
be —1.123 for v = 1 and —1.077 for u = 2. Both these values are well within the limits, viz.. 
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— 1.33 < < 0.79, obtained in model-independent estimates with the SN+CMB+BAO data 
in ref . [14] . On the other hand, the red-shift z = at which the best fit makes a transition 
from a value > — 1 to a value < —1 is found to be 0.2155 for u = 1 and 0.2450 for u = 2. 
However, stays well below zero even for z = 1.75, implying that DE is varying slowly with 
red-shift. The above values of Zc also agree fairly well with other independent studies [50]. The 
best fit DE density at the present epoch, p^^, is found to be equal to 0.7210 for z/ = 1 and 0.7184 
for u = 2. Remembering the dimensional re-scaling of the DE density, viz., p^/p^^, that 

we have performed earlier in Eq. (21), one may note that the p^^ shown in Fig. 1 is identical 
with the present DE density parameter Qqx = Pox/Poc (by virtue of the dimensional re-scaling). 
In other words, since the DE density p^ is effectively measured in units of the present critical 
density p^^, one has p^^ = fiox- It niay also be noted that the sum of the best fit Qq^ and the 
best fit flox is exactly equal to 1 (for both z/ = 1 and = 2), as it should be in accord with our 
prior assumption of the spatial flatness of the metric. This therefore proves the correctness of 
the x^-fitting of the model with the observational data. 




Figure 2. la, 2a and 3a contours in the parameter spaces A — B (for best fit J^om)) 
ilom — A (for best fit B), and Oom — B (for best fit ^4), are shown for the choices 
f = 1 (upper panels) and v = 2 (lower panels). The best fit points for both the choices 
are shown by the dots at the middle of all the la contours, whereas the cosmological 
constant, which corresponds to ^ = = 0, is shown by the dot that is found to lie on 
edge of the la A — B contour for both the choices (left panels, upper and lower). 

The la, 2a and 3a contour plots of (i) A versus B (with f2om fixed at its best fit value), 
(ii) versus A (with best fit B), and (iii) Qom versus B (with best fit A), are shown in Fig. 
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Figure 3. Left panels: extrapolations of w-j^,w and q (best fit with SN+CMB+BAO 
data), as functions of the scale factor a, to the range < a < 2, for the choices = 1 
(upper left) and v = 2 (lower left). Right panels: extrapolated variations of and 
(best fit with SN+CMB+BAO data) with the scale factor a, to the range < a < 2, for 
V = 1 (upper right) and v = 2 (lower right). The a = 1 line denotes the present epoch. 
The transition from a decelerating regime to an accelerating regime (i.e., the change of 
sign of q) takes place at a ~ 0.64. The dark energy density p^ is found to nearly follow 
the matter density p^ for a considerable period in the past until becoming dominant 
very recently, and increases slowly in the future. 

2, for the choices u = 1 (upper panels) and u = 2 (lower panels). The case A = B = 0, which 
resembles a cosmological constant DE, is found to be about la away from the best fit point in 
the A versus B contours (left panels), for both the choices. 

The upper and lower left panels of Fig. 3 depict the variations of the best fit DE EoS 
parameter w^{a), as well as the total EoS parameter w{a), Eq. (10), and the deceleration 
parameter q{a), Eq. (11), obtained as functions of the scale factor a (using the best fit values of 
the parameters flom, A, B) and extrapolated to the range < a < 2, for the choices z/ = 1 and 
V = 2 respectively. The range covers all of the past, i.e., right from the big bang (a = 0, 2; = 00) 
to the present (a = 1, 2; = 0), and a considerable part in the future, up to a = 2 (z = —1/2), i.e.. 
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when the present size of the universe gets doubled. Both and w are negative in the past and 
tend to become constant at a value close to each other and a little less than —1 in the future. 
The value of g, on the other hand, changes from positive to negative, i.e., the transition from 
deceleration to acceleration takes place at a = 0.640 {z = 0.562) for z/ = 1 and at a = 0.624 
{z = 0.603) for u = 2. In the future, q also remains negative and tends to be steady at a 
value close to w and . Thus the accelerated regime g < 0, as well as the 'super-acceleration' 
{w^ < —1), do not appear to be transient in the present model. 

The variations of the extrapolated best fit DE density and the matter density , with 
the scale factor a in the range < a < 2, are shown respectively for the choices u — 1 and v — 2, 
in upper and lower right panels of Fig. 3. For a considerable period in the past the DE density 
nearly follows the the track of the matter density, until exceeding the latter at scale factor 
a ~ 0.75, and dominant thereafter. In other words, decreases with a in a similar manner as 
p^ does in the early regimes, until at a recent epoch a ~ 0.75, when the DE begins to dominate. 
This behaviour, although not distinctly similar to that due to the tracker quintessence fields 
[18], may perhaps stand as a possible resolution to the coincidence problem [54]. One can, in 
fact, trace the similarity of the early universe profiles of p^ and p^ to the form of the chosen 
ansatz (30) for the field solutions and the resulting expression (33) for p^ . In the early epochs, 
i.e., for small values of a, the DE density p^ in Eq. (33) is dominated by the inverse power-law 
term ~ Aa~^ , similar to the matter density p^ = flQ^a'^. However, since ly < 3 and the best 
fit value of A is of the order of the best fit fiom, Px smaller than p^, and decreases less rapidly 
than the latter, for sufficiently smaller values of a. As a increases, the value of p^ eventually 
exceeds p^ due to the presence of the positive constant term (= A, given by Eq. (38)) in the 
expression (33) for p^. The Bin a term in Eq. (33), which is negative for a < 1 (i.e., past), is 
on the other hand, rather sub-dominant compared to A and does not play a very significant role 
either in the past or in near future. This is the reason why, the DE density p^ increases slowly 
and does not shoot up to very high values even at a scale factor as large as a = 2, giving rise to 
singularities in finite future. Admittedly, of course p^ ^ oo as a ^ oo due to the presence of 
the logarithmic term in p^. Thus, the extrapolations of the cosmological quantities using the 
best fit values of the model parameters, obtained in the red-shift range < 2; < 1.75, appear 
to hold for very distant past and future. 

In what follows, we integrate the expressions (31) numerically in the next section and use 
the values of fiom, A and B best fit with the data, so as to determine the variations of the scalar 
fields (f) and ^ with the scale factor a. We also reconstruct the potential V, given in Eq. (34), as 
a function of a, using these values of the parameters (ilom, ^, B), and finally, we work out the 
approximate analytic expressions for the functional variation of V with (p and ^, in the regimes 
a (distant past) and a < 1 (recent past). 
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5. Reconstruction of the scalEir potential 



Let us recall Eqs. (31), from which one can derive the following equations for the derivatives 
of the scalar fields (f) and ^ with respect to the scale factor a: 



Ho(f>'ia) = 



VAa-" + y/Aa-" + k 



(49) 



V2 a H{a) ' a H{a) 

where H{a) is as given by Eq. (37) or (39), in terms of the model parameters {Qom, A, B). 

Assuming the initial condition that = ^ = at a = 0, one may re-write the above 
equations in integral form as 



where 



Ho 0(a) 



I±{a) = 



[7+(a) - 7+(0)] , Ho e(a) 



Again, denoting = 
Ho 00 




da 



1 ± Wl + 



2ka'' 
~A 



[/_(a)-7_(0)] , (50) 



(51) 



Co at the present epoch (a 



Ho io 




/_(0)] . 



Prom Eqs. (50) and (52), one therefore finds 

0(a) /+(a)-/+(0) i{a) 



/_(i)-/-(o)- 



(52) 



(53) 



00 /+(!)- /+(0)' Co 
Now, in order to perform the integrations in I±, Eq. (51), one has to assign a particular 
value to the parameter k. Eliminating k from Eqs. (25) and (35) we find that the kinetic 
interaction Q satisfies the following cubic equation involving the parameters /?, 7 and B: 

- + 2{pB-2)Q + = 0. (54) 

Solving this equation and substituting the feasible root|| back in the relation (25), one finds k 
in terms of f3, 7 and B only. For B, we use its values best fit with the observational data, given 
in table 1 for the choices u = 1 and u = 2. For f3 and 7, wc recall that their values should 
be such that the condition 7"^ <^ <^ 1 is satisfied and the KIDQ Lagrangian (14) could 
emerge as an approximation to the two-field DBl Lagrangian (see the Appendix). Henceforth, 
assuming typically P — 0.01 and 7 = 10^, we find k — 0.012 for u — 1 and k — 0.004 for u — 2. 
Also, using the best fit values of the parameters ftom,^ and B, on which H depends, and 
performing numerically the integrations in I±, Eq. (51), we finally determine the variations 
of the normalized scalar fields 0/0o and C/Co with the scale factor a, ioT u — 1 and u — 2. 
Similarly, we also find how the quantity V — Vq, given by Eq. (34), varies with a, for the same 
choices of u. Such variations, extrapolated to the range < a < 2, are shown in Fig. 4. We 
observe that both 0/0o and C/Co increases with increasing a, however C/Co grows much faster 
than 0/00, whose variation gradually decreases with a (see the left panels of Fig 4). As such 



II Feasibility here implies that one should pick only that root Q which is real and greater than unity, so that 
the presumption of the positivity of k is ensured through the relation (25). 
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Figure 4. Evolutions of the normalized scalar fields (jj/ipo and ^/^o (left panels), and 
the scalar potential V minus its present value Vq (right panels), with the scale factor 
a, for the choices v = \ and = 2. Such evolutions, which are reconstructed using the 
best fit values of the model parameters Oom,^ and B, have been extrapolated to the 
range < a < 2. The vertical line a = 1 resembles the present epoch. 



the ratio ^/^q, which has been less than 0/0o in the past (i.e., a < 1), becomes greater than 
0/00 for a > 1 and grows to high values as we extrapolate it to far future. The potential V, on 
the other hand, decreases from a very high value (^ Vq, its present value) in the early epochs, 
reaches a minimum (< Vq) at some point a = in the past, and increases steadily thereafter. 
However, the entire profile ofV — Vo for both u = 1 and u = 2 (shown in the right panels of Fig 
4), is not symmetric about the minimum value Vm — Vq. In fact, the asymmetry is more when 
z/ = 2, rather than when u = 1. That is, the potential, after reaching its minimum, increases 
rather slowly for greater values of u. The values of and Vm (not shown in Fig. 4) could be 
calculated by extremizing the expression (34) for V: 



'{6-u)A' 
6~B 



1/^ 



V„. 



Vn 



{6-u)A 6 B 



+ 



In 



'{6-u)A' 
6~B 



(55) 
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For 1/ = 1, am = 0.6856 and Vm — Vq — —0.0948, whereas for u = 2, — 0.6558 and 
Vm — Vq — —0.0643. Therefore, the greater the value of u, the earher is the occurance of the 
minimum in the past, and the lesser is its value in magnitude. 

The overall variation of the potential V with the scale factor a could be explained as 
follows: In the very early epochs V — Vq, given by Eq. (34), is very large and positive due to 
the dominance of the positive inverse power-law term {6/u — l)Aa~''. With the increase of a, 
this term rapidly diminishes, and V — Vq becomes negative when the term 6B In a, which is 
negative for a < 1 (past), starts to dominate over the positive second term on the right hand 
side of Eq. (34). Eventually, V — Vo reaches the minimum, and then the potential V starts to 
increase as the term 6i?lna, though negative, gradually decreases in magnitude. V becomes 
equal to Vq at the present epoch (a = 1) and after that V — Vq increases with positive values as 
the logarithmic term becomes positive and increases with a. The asymmetry of the two sides 
of the minimum is obvious, because one is due to a power-law fall off and the other is due to a 
logarithmic increment. Also, for a bigger value of i> (here i> = 2), the power-law fall off is faster. 
The asymmetry is therefore more distinct, the minimum is attained earlier, and the minimum 
value Vm is smaller in magnitude. 

To reconstruct the potential F as a function of the fields and ^, we need to solve the Eqs. 
(49) (or, equivalently need to work out the integrals I±, Eq. (51)) analytically. However, this 
is very difficult because of the fairly complicated form of the normalized Hubble parameter H, 
given by Eq. (39). As an alternative, we resort to the following two regimes which are relevant 
for us: (i) a <^ 1 (early past) and (ii) a ~ 1 (recent past, present, and near future), and work 
out the approximate functional form of V{(j),^) in these regimes. 

(i) Ecirly Universe: 

For a ^ 1, the Hubble expansion is dominated by the inverse power-law terms in Eq. (37). 
As such, one can approximate: 

77-2/ \ , ^Om 

H^ia) fa + — ^ . (56) 

Now, from Eqs. (50) and (51), we have 

Hoim + e(a)]= V2A[I{a) - /(O)] , (57) 

where 

7(a) = -[J+(a)+7_(a)] - T . (58) 

Using the approximated form (56) of H, we get 

n«) ^ V^aI^ sinh-' ( a/ a^'-"^/' ] ^ 7(0) ^ , (59) 

whence 




Ho [0(a) + e(a)] « ^i^h"^ J a^'-^^'^ 1 . (60) 



3 — V 3 \ V i^^< 



Om 
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Inverting this relation and substituting in Eq. (34), one finally obtains 



+ 



6 



V 

12B 
3-u 



-I] A 



In < sinh 



3-z/ /3 



2 

3-u 



sinh-2-/(3--) 



3-z/ /3 



In 



Om 



6A 



(61) 
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Figure 5. Exact (solid line) and approximated (dashed line) variations of the quantity 
Hoi4' + with the scale factor a in the early epochs (a <^ 1), for the choices v = 1 
(upper left) and v = 2 (lower left). The approximation is found to hold up to a ~ 0.3 
for 1^ = 1, and up to a ~ 0.7 for v = 2. The right panels show the variation of the 
approximated V — Vq with i^o(</' + for the range of validity of the approximation, in 
the cases 1^ = 1 and v = 2 (upper right and lower right, respectively). 

The left panels of Fig. 5 show how the approximated form of the quantity Ho{(j) + as 
well as its exact form (obtained by working out the integral /, Eq. (58), numerically), vary with 
the scale factor a, for the choices u = 1 (upper left) and u = 2 (lower left). Whereas for v = l, 
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the approximation is found to be vahd only up to a ~ 0.3, it holds good till a ~ 0.7 for v = 2. 
Within the region of validity of the approximation, i?o(0 + increases almost linearly with a 
for = 1, whereas for z/ = 2, i/o(0 + O increases but gradually slows down as a increases. The 
variation of the approximated V — Vq as a function of the fields and ^, Eq. (61), is shown for 
V = 1 and i/ = 2 in the upper right and lower right panels of Fig. 5 respectively. Both these 
plots extend up to the range of a for which the approximation is valid in respective cases. The 
potential varies smoothly (i.e, without any discontinuity or multi valued- ness) with ifo(0 + O> 
as with the scale factor a (in Fig. 4). Also since i?o(0 + O increases monotonically with a, the 
nature of the V — Vq versus ifo(0 + plots in Fig. 5 is similar to the nature of the V — Vq 
versus a plots in Fig. 4 for smaller values of a. 

(ii) Recent Universe: 

Expanding the expression (39) for H in powers of (1 — a), for a 1 (i.e., close to the 
present epoch), and retaining only the terms linear in (1 — a) we have 



H\a) 



1 + hil 



where h = A - B + 



(62) 



for both y — 1 and v — 2. Now, from Eqs. (50), (51) and (52), one obtains^ 

Ho ma) - 0o) + iiia) - 6)] = V2A [7(a) - 7(1)] , (63) 

where 7(a) is the integral given by Eq. (58). One can evaluate 7(a) numerically for the choices 
V — 1 and V — 2 and find 77o [(0 — 0o) + (C ~ Co)], using the best values of the parameters 
ilom;^ a-nd B. However, to determine the functional form y(0, ^), we need to work out the 
integral 7(a) analytically. Let us separately consider the cases v — 1 and v — 2 diS follows: 

For V — V. The expression (34) for the potential V can be approximated as 

V{a) = Vo +{5A-6B){l-a). (64) 
On the other hand, the approximate analytic evaluation of the integral 7(a), Eq. (58), leads to 

1/2 

77o[(0(a)-0o) + (e(«)-6) 



2V2A 



1 + h 



1 + h 



-h 



- 1 



(65) 



Inverting this expression and substituting in the above equation (64), we get 



V{(p,0 ~ Vo +{5A-6B) 



1 - 



1 + h 



2V2A 



[(0 - 0o) + (e - eo)] +h 



(66) 



The variations of the exact and approximated forms of 77o [(0 — 0o) + (C — ■Co)] with the 
scale factor a, as well as the functional variation of the approximated [V^(0, C) ~ H], are shown 
in the upper panels (left and right respectively) of Fig. 6. The approximation is found to hold 
for a fairly large range 0.7 ^ a < 1.4. The approximated [\^(0, ^) — Vq], which has been plotted 



^ Note that Eq. (57) cannot be used now, because the approximation does not hold for a = 0. 
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Figure 6. Exact and approximate (solid and dashed lines) variations of 
Hq [{(j) — (po) + — ^q)] with the scale factor a(« 1), as well as the approximated 
functional variation [V{(l),0 ~ are shown for the choices 1^ = 1 (upper panels) 
and = 2 (lower panels). The approximation is found to be good in the range 
0-7 ^ a < lA for both v = 1 and 2. Within this range V — Vo has a minimum at 
Ho [{(p - 0o) + {(.- Co)] - -0.25 (for = 1) and ~ -0.15 (for iy = 2). 

for this range of validity, has a minimum (~ —0.1) at Hq [{(p — (Pq) + — ^q)] ~ —0.25. The 
overall profile of the approximated — Vq is similar to the exact variation oi V — Vq with a 
(shown in Fig. 4) in the range 0.7 < a < 1.4 where the approximation is found to be valid. 

For ly = 2: The approximation of the expression (34) for the potential V is given by 

V{a) = Vo + 2{2A-3B){l-a) , (67) 
and by evaluating the integral I (a) given by Eq. (58), approximately, one finds 



ifo[(0(a)-0o) + (e(a)-eo)] 



2A 
1 + h 




h 1 I ha 

1+h aV l+h 
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This expression cannot be inverted to get the scale factor a as a function of 
Hq — (po) + — Co)]) and hence obtain the expression for the analytic functional variation 
of the approximated V{(f),^). However, one may find the parametric plot of the approximated 
C) — Vo] versus Hq [{(f) — 0o) + ~ Co)] using the above equations (67) and (68). Such a 
plot is shown in the lower right panel of Fig. 6. [T^(0, ~ ^] varies in a similar way as for i/ — 1 
(see the upper right panel of Fig. 6), however the minimum value is higher and the minimum is 
reached at a value Hq [(0 — 0o) + (C ^ Co)] — —0.15, greater than that for z/ = 1. The lower left 
panel of Fig. 6) shows how the exact and approximate forms of Hq [{(f) — 0o) + (C ^ Co)] vary 
with the scale factor a. Similar to the case u = 1 (upper left), we find that the approximation 
remains valid in a fairly large region 0.7 < a < 1.4 as well for d = 2 (lower left). 

Finally, it should be mentioned here that although we have expressed the potential V 
approximately as a function of {(f> + C) for both a and a ~ 1 (see the Eqs. (61) and (66)), 
strictly speaking this cannot be true at all epochs. Indeed, if V is exactly a function of (0 + C), 
i.e., dV/d(f) — dV/d^, then the ansatz (30) becomes inconsistent with the scalar field equations 
of motion (19), given in sec. 3, and all our above results are invalid. The exact form of the 
potential V, which we are studying in a later work (in preparation) [54], should therefore be 
asymmetric in (f) and C (at least to their linear order) in order the model to be consistent. 

6. Conclusions 

We have thus explored the plausible crossing of the cosmological constant (A) barrier by the dark 
energy equation of state parameter in a fairly simple set-up of two canonical (quintessence- 
type) scalar fields with a mutual kinetic interaction. Such a crossing, which has been a big 
problem in many scalar field DE models, is shown to be realized with a specific form of the 
kinetic interaction and with the requirement that the dynamical part of total DE Hamiltonian 
is positive definite, so that the model is quantum mechanically consistent. Classical stability of 
the model is also guaranteed as the squares of the sound speeds corresponding to the adiabatic 
and entropy perturbation modes are positive definite. 

Under certain limiting conditions, the specific form of the kinetic interaction, which we 
study, can be shown to have originated from a higher dimensional two-scalar DBI action, that 
appears in the string theoretic scenario [44]. Such a kinetic interaction provides additional 
fiexibility in (apart from those provided by the usual kinetic terms of the scalar fields), so 
that the A-barrier (i.e., = —1) could be crossed at a particular epoch. 

Joint constraints on the parameters of the model by the SN+CMB+BAO data due to 
the SN Search Team [8], WMAP [14] and SDSS [15], show that has most likely crossed 
— 1 at a recent red-shift z = Zc (0.215 < Zc < 0.245), and its value w^^ at present is 
less than —1 (—1.123 < w^^ < —1.077). On the other hand, the transition from the 
decelerated phase of expansion of the universe to the accelerated phase takes place between 
0.562 < z < 0.603. All these results are fairly consistent with the model-independent estimates 
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with the SN+CMB+BAO data in ref. [14]. Additionally, we also observe that the dark energy 
density (best fit with the observational data) nearly follows the matter density (i.e., exhibits 
a similar fall-off with the scale factor a, as the latter) at early epochs, until exceeding it 
very recently. This apparently could provide a resolution to the coincidence problem (that is 
associated with the cosmological constant). Extrapolations to future epochs also show that the 
best fit dark energy density increases fairly slowly even at a fairly large scale factor, implying 
that singularities in finite future may plausibly be avoided in our model. 

The numerical reconstruction of functional forms of the scalar fields and the scalar 
potential, using the best fit values of the model parameters, shows smooth variations with 
the scale factor a, although analytically the exact form of the potential as function of the 
fields is very difficult to obtain. Working out therefore the approximate solutions for the scalar 
fields in the early universe and near the present epoch, we have obtained approximate analytic 
functional forms of the potential in terms of the scalar fields, in these regimes. Such analytic 
forms also exhibit the same smooth nature as the numerically reconstructed potential. 

Some interesting questions that arise in the context of the present model are in order: 

• Can one unify dark matter and dark energy in the general framework of a kinetically 
interacting double or multi-scalar theory, instead of treating them separately as in this 
model? 

• Can we generically determine for kinetically interacting double or multi-quintessence 
model, the exact form of the scalar potential, which could lead to the A-barrier crossing as 
well as the dark matter tracking (by DE)? If so, then could it be ascertained whether such 
a potential belongs the class of tracking or scahng potentials that arise in generic k-essence 
theories? 

• Can we have in the context of a kinetically interacting double or multi-quintessence model, 
the assisted accelerated solutions, which have been shown to exist generically for the multi- 
field k-essence models admitting scaling solutions [55]? 

• Can we ascertain the status of the future singularities, if any, in the context of kinetically 
interacting double-quintessence models, generically, i.e., by not just resorting to a 
particular ansatz to solve for the field equations? Or, can we generically ascertain whether 
or not the cosmic super-acceleration {w-^ < — 1 regime) is always eternal (as in the present 
model) for kinetically interacting two-field quintessence? 

Works addressing some of these questions are in progress [54], which we hope to report soon. 
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Appendix: Kinetic interaction from a two-field DBI perspective 

The Dirac-Born-Infeld (DBI) multi-scalar Lagrangian is particularly important in view of the 
notion acquired from string theory that our observable four dimensional world may be looked 
upon as being a warped D3-brane embedded in a higher dimensional (bulk) space-time [56]. 
The general expression for such a Lagrangian in an effective four dimensional theory is given 
by [41, 42]: 

P = 7 - - U{4)') , (69) 

where (f)\l = 1,2,...) are a set of scalar fields, which from the bulk point of view, correspond 
to the coordinates of the brane in the extra dimensions, and U (0^) is the multi-field potential 
due to the interaction of the brane with bulk fields or with fields on other branes. T> is the 
determinant of induced metric on the brane, given by 

V = det {5; + Gu d^^'d'cp') , (70) 

where Qu is the field space metric, which is proportional to the extra dimensional metric living 
in the bulk, and 7 is a coupling parameter, that appears by virtue of the warping of the D3- 
brane in the bulk. Both and 7 could generally be functions of the fields 0^, and in a 
homogeneous FRW background the above expression for V takes the form [41, 42]: 

V ^ 1 - 2r'gijx'', x^'^^. (71) 

Now, for a configuration of two fields 0^ := {4>,i}, assuming 7 and Qu to be constants 
(for all I, J), we can write 

2^ = 1 - 7"' (^^11 0' + G22 + ^Gi2 0e) ■ (72) 

Rescaling and ^ such that both the metric components Qu and Q22 are effectively set to unity, 
we get 

T> — 1 — ^ ^ — —0^ , where /? = = constant . (73) 
72 7 

Under the assumption: 

7 > ^12 > 1 /3 < 1 , but 7-^ < /3 , (74) 

one can write"*" 



Substituting this in Eq. (69), and neglecting the O {P/'j) and higher order terms, we finally 
obtain 



P = - 7 yi - - (76) 

+ Note that in this paper we have typically set 7 = 10* and with Q12 = 10^, it implies (3 = 4^12/7 = 10~^ <C 1 
(but 7~^), as per the assumption (74) [see sec. 5]. 
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This is the same as the KIDQ Lagrangian (14) considered in sec. 3; V{(j),^) = f/(0, ^) — 7 being 
the effective (shifted) two-scalar potentiaL 

Stability criterion: In general multiple scalar field models, the cosmological perturbations 
of scalar type are divided into: (i) the adiabatic (instantaneous) modes, which are fluctuations 
along the field space trajectory, and (ii) the entropy modes, which are orthogonal to the field 
space trajectory [57]. The squares of the speeds of propagation of both these modes should be 
positive definite in order that the underlying model is cosmologically stable. 
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Figure 7. Plots of the KIDQ sound speed Cg (assumed to be isotropic) as a function 
of the scale factor a, for the choices v = 1 and v = 2. Except at very early epochs 
(a <C 1), the sound speed is found be extremely close to unity (i.e., the speed of light), 
for both the choices. 

For multi-scalar DBI models [41, 42], the adiabatic and entropy modes are shown to be 
isotropic, i.e., they propagate with the same (sound) speed equal to \/V. Since the KIDQ 
Lagrangian is an approximation of the two-scalar DBI Lagrangian under the assumption (74), 
one may intuitively consider this isotropy of the perturbation modes exists in the case of KIDQ 
as well, and that the effective sound speed for KIDQ to be Cg = \/V, given by Eq. (75). 
Obviously, as shown in Fig. 7, this sound speed (and as such its square) is positive definite and 
very close to unity (speed of light) for both v = 1 and = 2, since the 7"^ and (3 terms in Eq. 
(75) are very small compared to unity. 

More generally of course, if we treat KIDQ as an exact model (i.e., not an approximation 
of DBI), then a rigorous analysis following the general formalism worked out in [42], shows 
that the sound speeds corresponding to the adiabatic and entropy modes do actually differ. 
However, the difference is very slight for the values of the parameters /3, 7 used in this paper, 
and the sound speeds are still positive definite and close to unity [54]. 
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